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We investigate the reach of the LHC Run 2 and that of a future circular hadron collider with up 
to 100 TeV centre of mass energy for the exploration of potential Dark Matter sectors. These 
dark sectors are conveniently and broadly described by simplified models. The simplified models 
we consider provide microscopic descriptions of interactions between the Standard Model partons 
and the dark sector particles mediated by the four basic types of messenger fields: scalar, pseudo¬ 
scalar, vector or axial-vector. Our analysis extends and updates the previously available results for 
the LHC at 8 and 14 TeV to 100 TeV for models with all four messenger types. We revisit and 
improve the analysis at 14 TeV, by studying a variety of analysis techniques, concluding that the 
most discriminating variables correspond to the missing transverse energy and the azimuthal angle 
between jets in the final state. Going to 100 TeV, the limits on simplified models of Dark Matter 
are enhanced significantly, in particular for heavier mediators and dark sector particles, for which 
the available phase space at the LHC is restricted. The possibility of a 100 TeV collider provides an 
unprecedented coverage of the dark sector basic parameters and a unique opportunity to pin down 
the particle nature of Dark Matter and its interactions with the Standard Model. 


I. INTRODUCTION 

The data collected by the Planck mission [1] confirms that, based on the standard model of cosmology, dark matter 
constitutes nearly 85% of the total matter content in the universe. With a natural assumption that all matter in 
the universe, dark and visible, is fundamental, dark matter should be described by a microscopic particle theory^. 
Hence the quest to establish the identity of dark matter, and its fundamental interactions, amounts to one of the 
most important goals in particle physics. 

The observational evidence for dark matter (DM) was established from gravitational effects on visible matter. 
However, the Standard Model (SM) of particle physics does not contain any viable DM candidates. In this way dark 
matter provides us with arguably the strongest experimental evidence for the existence of physics beyond the Standard 
Model. The observation of non-gravitational interactions of DM with visible matter could be crucial in discovering 
extensions of known fundamental theories. Since the experimental evidence for the existence of DM is currently only 
gravitational, at present we have no definitive microscopic understanding of dark matter. For these reasons the search 
for DM has escalated in recent years with direct detection, indirect detection, and most recently high energy collider 
searches combining to restrict the range of parameter space for dark matter models. 

Presently the frontier of high-energy physics corresponds to the newly upgraded LHC, which is currently operating 
at ^/s = 13 TeV, with a planned increase to y/s = 14 TeV in the near future. However, it is certainly plausible that 
the scale of new physics is out of range of the LHC’s reach. Broadly speaking, analyses of Run 1 data are able to 
constrain new physics approximately up to 2 TeV. The higher operating energy at Run 2 will certainly increase these 
bounds, but not by more than an order of magnitude. Therefore if new physics arises first at scales greater than those 
probed by the LHC, it may not be possible to infer its existence with the current generation of machines. For this 
reason, the community is beginning to consider the potential for the successor to the LHC [3j. The next machine, will 
have to significantly extend the reach of the LHC, and as such, requires a large increase in operating energy. Current 
proposals set 100 TeV as the target center of mass energy. Since a hadron machine of this energy mandates a circular 
design the putative machine is usually referred to as the Future Circular Collider (FCC). Such a machine would be a 
formidable tool in the quest to discover new physics. 

The main motivation of this paper is to investigate the reach of the FCC with a 100 TeV centre of mass energy for 
dark matter searches and ideally potential discoveries. In parallel with our earlier work [3] which analyzed the limits of 


^ for a review see e.g. [2] 
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the LHC at 8 TeV and its DM discovery potential at 14 TeV, we aim for a general mo del-independent characterisation 
of dark sectors with as few free parameters as possible, but without resorting to an effective field theory description. 
To achieve this we employ the simplified model approach which provides a microscopic QFT description of a minimal 
set of interactions between the Standard Model partons and the dark sector particles. These interactions are mediated 
by a complete set of four basic types of messenger fields, i.e. the scalar, pseudo-scalar, vector and axial-vector. One 
naturally expects that the FCC should perform significantly better than the LHC for heavy mediators and heavy 
dark matter particles, the primary aim of this work is to seek to quantify this improvement. Our work is another 
step in the emerging program of DM studies at future colliders in the 100 TeV range, Related studies using 

simplified models for constraining dark sectors at the LHC include Refs. [SHU], and we also refer the reader to the 
recent summaries [171118j and references therein. 

In DM searches at hadron colliders, the putative dark particles are pair-produced in collisions of the visible sector 
particles - the Standard Model quarks and gluons. In the set-up we study here, there are no direct interactions 
between the SM sector and the dark matter particles. Instead these interactions are mediated by an intermediate 
degree of freedom - the mediator field. In general, one can expect four types of mediators, scalar S, pseudo-scalar P, 
vector Z' or axial-vector Z". The corresponding four classes of simplified models describing elementary interactions 
of these four mediators with the SM quarks and with the dark sector fermions x ^'I'e 
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The coupling constant ^dm characterizes the interactions of the messengers with the dark sector particles, which for 
simplicity we take to be Dirac fermions x, X, tti® case of scalar DM particles is a straightforward extension of these 
results. 

The coupling constants linking the messengers to the SM quarks are collectively described by 


scalar & pseudo — scalar messengers : 
vector & axial — vector messengers : 
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For scalar and pseudo-scalar messengers the couplings to quarks are taken to be proportional to the corresponding 
Higgs Yukawa couplings, yq as in models with minimal flavour violation |I9j . and we keep the scaling gq flavour- 
universal for all quarks. For axial and vector mediators gsM is a gauge coupling in the dark sector which we also take 
to be flavour universal. The coupling parameters which we can vary are thus ^dm plus either gq or ^sm, the latter 
choice depending on the messengers.^ 

In general, the simplified model description of the dark sector is characterised by five parameters: the mediator 
mass tomed, the mediator width Fmed, the dark particle mass muM, and the mediator-SM and the mediator-Dark 
sector couplings, 5sm, guM- Out of these, the mediator width Fmed, does not appear explicitly in the simplified 
model Lagrangians Q-Q and should be specified separately. Fmed accounts for the allowed decay modes of a 
given mediator particle into other particles from the visible and the dark sector. In a complete theory, Fmed can be 
computed from its Lagrangian, but in a simplified model we can instead determine only the so-called minimal width 
Fmed. min, i-c. the mediator width computed using the mediator interactions with the SM quarks and the x, X DM 
particles defined in Eqs. 0-0- Importantly Fmed, min does not take into account the possibility of the mediator to 
decay into e.g. other particles of the dark sector, beyond x, X, which would increase the value of Fdm- In Ref. |3] 
we have investigated the role of F med as an independent parameter in the simplified models characterisation of dark 
sectors by using a simple grid for Fdm = {1, 2, 5,10} x Fmed, min- 


^ In Ref. [4] we have parametrised ^dm for (pseudo-)scalar messengers as poM = 5 x™dm/'d to look symmetric w.r.t. 0 , and have 
treated as a free parameter. Here we do not impose this requirement and leave guM the free parameter. 



3 


Here we will not repeat this analysis, referring to [3] for a more general discussion of Fmed; we will instead adopt 
a reduced simplified description where the width is set to its minimal computed value Fmed. min which amounts to 
larger signal cross-sections (we will also check that Fmed, min < ?^med/2). For our simplified models we have 


Nf 

Fmed, min — ^XX ^ ^ -^c F 
i=l 


(7) 


where F^^ is the mediator decay rate into two DM fermions, and the sum is over the SM quark flavours. Depending 
on the mediator mass, decays to top quarks may or may not be open i.e. uimeu should be > 2mt for an open decay. 
The partial decay widths of vector. Axial-vector, scalar and pseudo-scalar mediators into fermions are given by. 
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where mj denotes masses of either SM quarks q or DM fermions x the coupling constant gf denotes either gsM 
or 5 dm- 




FIG. 1: Representative Feynman diagrams for gluon and quark induced mono-jet plus MET processes. The mediator X can 
be a scalar, pseudo-scalar, vector or axial-vector particle. The gluon fusion process involves the heavy quark loop which we 
compute in the microscopic theory, while the quark-anti-quark annihilation is a tree-level process at leading order. 


In this paper we will focus on jets plus missing energy searches, generalizing our earlier results [3] from the 8-14 
TeV to 100 TeV colliders. The analysis in [4] was based on a mono-jet plus missing transverse energy (MET, or ) 
signature - a popular choice used in searches for new physics including supersymmetry, extra dimensions and dark 
matter at the Tevatron and the LHC [IDHSlj. Here we will update the experimental analysis techniques to take into 
account both the leading and trailing jets in the final state to present a more realistic idea of potential limits. Hence 
we will update our LHC 14 results j4] accordingly to provide a fair benchmark. 

Depending on the choice for the mediator field different production mechanisms will contribute. For vectors and 
axial-vectors the dominant mechanism is the quark-antiquark annihilation at tree-level. For scalars and pseudo-scalars 
on the other hand, the loop-level gluon fusion processes are more relevant. The representative Feynman diagrams for 
both channels are shown in Fig. [l] In comparing DM collider searches with direct and indirect detection experiments 
it is important to keep in mind that our collider processes and limits continue to be applicable for discovery of any dark 
sector particles escaping the detector. Hence dark particles produced at colliders do not have to be the cosmologically 
stable dark matter. 

Finally we would like to comment on the possible origin and the UV consistency of the simplified models 0-0- 
The scalar and pseudo-scalar messenger fields in our simplified models 0-0 are singlets under the Standard Model. 
How can this be reconciled with the fact that they are supposed to be Higgs-like, with the Higgs being an SU{2)e 
doublet? In fact, the simplified models 0-0 can arise from two types of the more fundamental theories. The 
simplest theories of the first type are the two-Higgs-doublet models [33]. In this case the mediators would originate 
from the second Higgs doublet. The other type of models giving rise to our simplified models are even simpler in the 
sense that scalar mediators (and the dark sector particles they are coupled to) can be genuinely neutral under the 
SM but mix with the neutral component of the Higgs. Following the Higgs discovery there is a renewed interest in 
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FIG. 2: Reconstruction efficiency for signal Csig and background eskg as functions of individual kinematic variables at 14 TeV 
with 1 ab“^ luminosity for scalar mediator models using single variables as well as the combinations of variables, as indicated. 
The dashed lines depict single jet signal MC, while the solid lines correspond to merged samples with np to two jets. 


the literature in Higgs portal models [31H37] where the scalar mediators apartonsre SM-singlets but the SM Higgs h 
interacts with them via the interaction, AhplHpl^p. When and H both develop a VEV, mass mixing occurs, and 
after transforming into the mass eigenstate basis, one finds two scalar resonances hi and ^ 2 , both of which interact 
with the Standard Model and the dark sector, with the hi state identified with the SM Higgs and /12 being the scalar 
mediator. These models provide a direct connection of the dark sector with Higgs physics and can link the origin of the 
electroweak and the DM scales [HHl - HT] . The simplified dark sector models with vector and axial-vector mediators in 
Eqs. (|^-Q can also be derived from appropriate first-principles theories. Since the mediators are spin-one particles, 
these UV models would necessarily require the mediators to be gauge fields and the DM to be charged under these 
gauge transformations. A classification of anomaly-free extensions of the Standard Model Abelian 17(1)' factor was 
given in |42j and can be used for constructing an example of a consistent gauge-invariant vector and axial theories of 
the type Q. 

This paper proceeds as follows, in section [H] we discuss the theoretical and experimental setup we used for our 
analysis. The results of our analysis are presented in section m and our conclusions are presented in section [TV] 


II. THEORETICAL AND EXPERIMENTAL SETUP 
A. Collider design and Monte Carlo generation 

At this moment, the exact nature of the FCC is unclear. This level of uncertainty includes fundamental issues, like 
the operating center of mass energy of the machine, and technical details including, for instance, the capabilities of 
future detectors. This makes robust predictions impossible at the present time, instead it is more interesting to study 
what is potentially feasible at future colliders, given a modest (and conservative) set of assumptions. 

We begin by discussing the straightforward benefits associated with higher energy machines in relation to DM 
searches. Clearly with larger center of mass energies, heavier mediators can be produced on-shell, allowing for 
significant enhancements for heavy mediators compared to a 14 TeV analysis. Secondly there are effects due to the 
parton distribution functions (PDFs). A typical 14 TeV DM analysis naturally focusses on the region of high missing 
transverse energy, which corresponds to larger values of the partonic center of mass energy (s). This region of phase 
space is suppressed by the PDFs, and therefore additionally damps the differential cross section, beyond the simple 
phase space/matrix element falloff. However, at higher energies the same region of phase space now corresponds to 
much smaller s/s values, and as a result the PDF suppression is no longer as severe. Therefore, even for mass scales 
which are currently probed at the LHC, we expect a higher energy machine to be able to improve upon the bounds. 

A major concern for the simulation of analyses at an FCC, is the modelling of the theoretical event generation. It is 
clear that a huge amount of work will be needed over the coming decades to ensure that the Monte Carlo simulations 
provide a reasonable description of signal and background processes [l^. We do not attempt this huge feat in this 
paper, but instead list the concerns regarding existing tools, and their potential impact. 
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FIG. 3: Kinematic distributions for signal scalar and mediator models and the SM backgrounds at 14 TeV assuming 1 ab”'^ 
of integrated luminosity. We show four kinematic variables: pr of the leading jet, px of the trailing jet, the azimuthal angle 
between the jets cos{A(j>jj), and the missing energy . Ratios of {S + B)jB are shown for each observable. The red bands 
indicate the uncertainties on the background distributions. The accordingly color-coded numbers for sig(JJ) and sig(J) give 
the statistical significance to disfavour the presence of the signal using the CLs method. 


Firstly, there is the obvious issue regarding tuning and PDF fits, which have been undertaken at much lower energy 
than those accessible at FCC’s. Of particular importance are the extraction of the gluon PDF’s, which will dominate 
the initial state configurations. Secondly, there is the accuracy of the perturbative component of the simulations. This 
is particularly worrisome. At FCC energies emission of additional radiation will result in copious jet-production around 
the Electroweak scale. This will require delicate handling with respect to matching and merging of parton shower 
and matrix element emissions. Since matching prescriptions typically require scales which separate emissions into 
the two categories, one can easily imagine existing tools for LHC physics are not optimal for future FCC predictions. 
An additional concern relates to the simulation of Electroweak bosons, at 100 TeV, the mass of the W, Z and H 
bosons become small scales, and Sudakov logarithms associated with their emission from partons becomes relevant. 
The resolution of these issues will take many years of research and improvements to existing tools. However given the 
likely timescale of construction, and the rapid improvement in theoretical tools, none of the above issues should be 
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FIG. 4: Kinematic distributions for signal vector and axial-vector models and the SM backgrounds at 14 TeV assuming 1 ab”'^ 
of integrated luminosity. We show four kinematic variables: pr of the leading jet, px of the trailing jet, the azimuthal angle 
between the jets cos{A(j>jj), and the missing energy Jfix ■ Ratios of {S + B)jB are shown for each observable. The red bands 
indicate the uncertainties on the background distributions. The accordingly color-coded numbers for sig(JJ) and sig(J) give 
the statistical significance to disfavour the presence of the signal using the CLs method. 


regarded as significantly likely to negatively affect the physics program at the FCC. They should certainly be kept in 
mind and the resulting theoretical predictions used in this paper should be interpreted as having a large uncertainty. 

In order to attempt to simulate events at the FCC it is therefore critical to include high energy jet radiation in the 
best manner possible. This can be achieved by ensuring that matched samples are used. Backgrounds are generated at 
next to leading order for 0,1,2 jets merged using MadGraph-aMC@NLO [H], with the exception of FF-1-jets which does 
not have the second jet merged. Diboson and single top processes are all included. For the signal we use MadGraph 
for the Vector/Axial simplified models and a combination of MCFM |45l|46] and VBFNLO [47] for the production 
of Scalar/Pseudoscalar mediators in association with one and two-jets. The output THE events are then merged 
using the GKKW-L interface of Pythia 8 |1S|. NNPDF3.0 |15| PDF’s are used for the generation of all Monte-Garlo 
samples. In order to generate the samples consistently the mediator is assumed to be on-shell for the 1 and 2 jet 
merged sample. The one-jet sample is generated with the full width and thus includes off-shell effects. Since the 
putative mediators scrutinized at the FCC will include heavier masses, off-shell effects will indeed play an important 
role in the simulation. However, given the large center of mass energy, including at least the two-jet matrix element 
(and ideally in the future 3-jet and beyond) is likely to be a more important requirement than including off-shell 
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FIG. 5: Kinematic distributions for the two most discriminating variables and cos{A(j)jj) at = 100 TeV. 


effects^. In the following section we therefore present results both for the one-jet sample (with the full off-shell effects) 
and the two-jet merged sample (which includes hard matrix element jets, but requires on-shell mediators). We note 
that that this is not an issue for the vector/axial mediators. In all cases we calculate the minimum width which 
corresponds to all open fermionic decay channels (see Eqs. 0 -® for definitions). 

In addition to the theoretical issues discussed above, we have limited knowledge of future detectors. However, here 
a conservative assumption can easily be made. We simply assume in this paper, that the detector at the FCC is a 
copy of the CMS detector. One extension to this detector is added, the lepton acceptance for this upgraded detector 
is assumed to extend up to I 77 I < 4.0 and \r]\ < 5.5 for the 14 TeV and 100 TeV detector respectively. This design is 
consistent with the plans for the upgraded CMS and ATLAS detectors [snug. The choice of detector allows us to 
utilize the extensive studies of the detector performance. In reality we can expect improvements in technology which 
can significantly improve on existing detectors. 

To emulate the effects of a detector, the MET is smeared as a function of the mediator px [53]. The jets are smeared 
as a function of rj and px [M]. Effects from pileup are taken to account, such that the respected observed resolution 
follows closely the expected conditions for high luminosity running at the LHC. 


^ We note that very recently ref [50] presented an implementation of the signal model which can incorporate one- and two-jet merged 
samples and off-shell effects. 


























































B. Kinematic distributions for signal and background 


In light of using multi-leg Monte Carlos and recent developments in shape based analysis for Dark Matter projec¬ 
tions [55] , the analysis selection was re-optimized from a previous cut and count analysis by considering a shape based 
analysis in both single and pairwise combinations of variables that take into account both the leading and trailing 
jets. To hnd sensitive variables to separate signal from backgrounds we focused on , rri^, the razor variables, 
the scalar sum of all jets above 30 GeV (/it), and the angle between the two leading jets jSSHSTj. Figure]^ 

shows the performance of these variables when compared with an amalgam of signal Monte Carlos. For these signm 
models, a dark matter mass of 50 GeV is taken (ensuring on-shell decays for both). The range of samples consist of 
those with different mediator masses ranging from ITeV to 2 TeV in 100 GeV intervals. The samples are combined by 
cross-section weights. We recall that the 1 jet Monte Garlo is generated with full on and off shell width effects, whereas 
the l-jet/2-jet merged sample is generated explicitly on-shell. When considering the single variable performance, we 
hnd is the most discriminating variable. As a result, combining these variables into pairwise combinations using 
a boosted decision tree we hnd that the most sensitive combination is and between the two leading jets. 

This can be readily seen in the left plot in Fig. where the A(j)jj+^x combination is shown as the solid line in blue. 
In the lower-left quadrant of the plot the reconstruction efficiency of the signal Cgig signihcantly exceeds that of the 
background, for example for Csig ~ 0.5 we have ebkg — 0.1 on the blue contour in Fig. A comparison with higher 
dimensional sets of variables in the right panel shows even further improvement, particularly when the second jet 
merged MG is used. However, for this further optimization to manifest itself in the analysis requires precise knowl¬ 
edge of the kinematic discrimination of in the far tails of the distribution. Additionally we note that the larger 
differences between the one-jet and two-jet signal MG, can be indicative of greater sensitivity to theoretical modeling. 
Thus, we do not extend these studies beyond their current projections. Our hnal extractions are then performed by 
fitting the distribution in two bins of A4)jj {A(pjj <1.1 and Acfijj > 1.1). The binning is optimized such that the 
the full spectrum is covered with MG samples. Additionally a selection of min(A0^^ > 0.5, where i runs over all 
jets is applied to minimize the impact of QCD multi-jet events where one or several of the jets are miss-measured. As 
a result of this selection being applied, this background is neglected throughout the course of these studies. Further 
improvements could be obtained by using hner bins for the tails of the distribution^. 

The signal extraction is performed with a full shape analysis in which the shape for the dominant backgrounds are 
taken from control regions. For the Z ^ vv background, the Z —> control region is used as a background. For 

the W —>■ tv, top and diboson backgrounds, we use the single lepton control region. For each of these control regions 
the full statistical uncertainty on the shape is propagated per bin on each of the backgrounds with an additional 
one percent uncertainty uncorrelated per bin to account for additional modelling uncertainties. For all but the tail 
bins of the shape uncertainties on the l/trp spectrum are roughly 1% with the dominant uncertainty resulting from 
the additional one percent modelling uncertainty. Finally, the signal is profiled using the standard limit extraction 
(GLg) [^[^. Additional nuisances are placed on the background normalization for lepton efficiencies and luminosity. 
These are constrained to very small values due to the large dataset in the signal region and do not affect the limit 
sensitivity. 

Kinematic distributions for our simplified models of dark sectors alongside the main SM backgrounds are shown in 
Figs. [^ [^ and The distributions are shown as functions of four kinematic variables, pt of the leading jet, pt of the 
trailing jet, the azimuthal angle between the jets cos{A(f>jj), and missing energy l^rp . 

In Figs. 1^ and 1^ we show shapes of kinematic distributions corresponding to the (pseudo-)scalar and (axial-)vector 
case respectively, with the signals and backgrounds computed at 14 TeV centre-of-mass enercy and assuming 1 ab“^ 
integrated luminosity. The only event selection cuts imposed for the distributions in Figs. m are > 200 GeV and 
min(A(/)^^ > 0.5, where i runs over all jets in each event. 

We hnd that for the event selection cuts imposed, the distributions for samples using multi-jet merging are similar 
to the unmerged samples. In the gluon-induced scalar and pseudoscalar cases the unmerged samples result in a harder 
Pt distribution, particularly for the trailing jet. On the other hand, the cos{A(j)jj) distribution changes more rapidly 
when including multi-jet merging. In Figure we show the kinematic distributions for the two most discriminating 
variables, which, as was pointed out above, are the missing transverse energy ^t the azimuthal angle between 
the jets A(j>jj at ^/s = 100 TeV. The two plots on the left show the vector and axial-vector case while the two plots 
on the right show the scalar and pseudoscalar signals. 
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While we will not include this approach in our analysis we point out that this could be a viable way to optimize the sensitivity of the 
search. 
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C. Direct and Indirect Detection Limits 


In Sec. Ill we will compare collider limits with limits from direct and indirect detection experiments. For convenience 
we recall the formulas for direct and indirect detection limits, in terms of the simplified models of Eq. 0-0- 
We use 


and 
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’’’ ’^MED 


3 

”^MED 


( 10 ) 


( 11 ) 


with a ~ 0.43 [501IST] and the reduced mass p = ■mDMWp/(TODM F^p), for the cross section of a dark matter particle 
scattering spin-independently (vector mediator) or spin-dependently (axial-vector mediator) from a proton. 

The cross section for a Dark Matter particle scattering from a nuclei via a scalar mediator of Eq. Q is given by 
[62H64] 


cr 
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TT 


E 

q—u,d,s 


■p 


f gDMgqUq \ Wp / gPMgqUq 

V ”^MED / 27 q=c,b,t ^ ”^MED 


( 12 ) 


where = 0.019, = 0.045, /f = 0.043 and /tg — 1 — J2q=u d s fq ESI ES] and mp is the proton mass. 

When comparing the expected sensitivity for the LHC and FCC for DM searches to those of Direct Detection 
it is interesting to compare to expected impact of the neutrino wall [iniET]. The neutrino wall is a result of the 
background of cosmic neutrinos. We take their interaction cross section to be indicative for the ultimate reach of 
DD experiments. Therefore in the following section we will present the equivalent limit which may be obtained from 
DD experiments with the currently hypothesized values of the neutrino wall uniET]. For a pseudo-scalar mediator, 
taking existing limits into account [681169j , indirect detection experiments can result in stronger limits than direct 
detection experiments [701 El]. For the simplified model of Eq. 0- using the velocity-averaged DM annihilation cross 
section into bb, 
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which allows us to derive a limit on the parameters in the bb channel 


(13) 


III. RESULTS 

Results are obtained scanning over a spectrum of signal models at 14 TeV and 100 TeV. A predicted luminosity of 
1 ab“^ is used for both analyses, so the sensitivity can be compared directly. We note that this amount of integrated 
luminosity is a rather modest amount compared to what is likely to be collected at a future collider, whereas for the 
LHC it corresponds to a significant fraction of the data set that will be obtained. As a result the results we present 
here should be interpreted as those which can be obtained over the lifetime of the LHC, and for a shorter run with 
the FCC. 

We begin by studying constraints on total cross sections which can be obtained using our analysis. Figure [6|presents 
the total cross section which the analysis excludes for each of the four mediator types defined in Eqs. 0^Q. We 
define our cross sections by setting ^dm = ffSM = 1 and select the mediator mass as indicated in the legend of each 
figure respectively. As an illustrative example we have chosen a relatively small characteristic value of 100 GeV, 
although the results obtained for other kinematically accessible values of dark matter mass were found to be similar. 
The kinematics of the process are then completely specified once the couplings (7 dm and gsM are set, since this fixes 
the minimal width of the mediator [4]. The excluded cross section is then related to the predicted cross section as 
follows, 


cr = p cr{guM = 1,5SM = IjTOmed), (14) 

With the kinematics of the model fixed we set a limit on p, defined above using the CLg-method, again assuming 
1 ab~^ of data. Values with p < 1 indicate the excluded couplings and width are smaller than the tested model, and 
the point is then excluded. 
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FIG. 6: Cross section exclusion limits as a function of mediator mass for a fixed dark matter mass at a given coupling. We 
show results for vector (upper left panel), axial-vector (upper right panel), scalar (lower left panel) and pseudoscalar (lower 
right panel) 


In Fig. 1^ we also distinguish between the mono-jet (shown in green) and the multi-jet-based analyses (shown in 
yellow). It can be seen that the new multi-jet-based analysis is more powerful and provides a considerable improvement 
at 14 and at 100 TeV. It readily follows from Fig.j^that the new multi-leg analysis provides a considerable improvement 
at 14 TeV relative to the results of our earlier work [3]. At 100 TeV the exclusion limits in Fig. are completely 
new, and we also point out that the /r-factor remains ^ 1 well into the TeV range of mediator masses even for scalar 
and pseudo-scalar mediators. The impact of including the additional matrix elements in the signal simulation is, as 
predicted, much greater for the FCC, which allows for copious production of light degrees of freedom. For the case of 
scalar and pseudo-scalar mediators at 14 TeV there is a cross-over for mediators heavier than ~ 1 TeV, which is absent 
at 100 TeV. This corresponds to exactly the regions of phase space in which the off-shell effects dominate. The one-jet 
sample has access to the significant cross section which arises from the tails of the Breit-Wigner distribution, whereas 
the multileg sample does not. This region therefore has large theory errors using the multi-leg sample. However, 
we note that the region of phase space for which the multi-leg sample breaks down is far from the values of p, = 1, 
so this region of phase space is of limited importance in regards to setting limits on model parameters. Finally we 
note that Fig. also includes cross sections for interesting SM predictions which the FCC and Run II of the LHC 
will investigate. We present the cross sections for ttH and HH and show their relative size compared to our DM 
predictions. 
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FIG. 7: Exclusion limits on mediator mass versus dark matter mass for given couplings based on multijet measurements at 
14 TeV LHC. We show results for vector (upper left panel), axial-vector (upper right panel), scalar (lower left panel) and 
pseudo-scalar (lower right panel) mediators. Plots for vector and axial-vector mediators also show the regions inside which 
Fmbd, min < uiMBD SO that the particle physics interpretation of (axial)-vector mediators applies. Scalar and pseudo-scalar 
mediators widths we use are sufficiently narrow. 


Figurej^shows the 14 TeV LHC exclusion limits for all four mediator models interpreted on the (mMED, wdm) plane, 
and the color-coding indicates the value of the mediator coupling, assuming = <7sm- To provide a benchmark 
for the collider reach for searches of DM, we compare their projected exclusion contours with those relevant to the 
present and future direct and indirect detection experiments. Specifically, for the cases of scalar, vector and axial- 
vector mediators we show the neutrino wall which is expected to correspond to the ultimate reach of the direct 
detection (DD) experiments. For the pseudo-scalar mediators one should compare with the indirect detection (ID) (to 
avoid the non-relativistic velocity suppression affecting pseudo-scalar cross-sections in DD experiments). Hence for 
the pseudo-scalar mediators we show the indirect detection limits using FERMI-LAT data [SS]. Figure [^generalizes 
the results to the FCC case with 100 TeV energy. We find large improvements in reach for the mediator mass ttimed 
as well as the dark matter mass todm- Increasing the center-of-mass energy from 14 to 100 TeV allows us to probe 
4-5 times larger mediator and dark matter masses, irrespective of the mediator’s spin. It is interesting to note that 
the limits obtained for the neutrino wall for vector, axial and scalar mediators included completely those of the LHC 
operating at 14 TeV. That is to say that, given enough time, data and understanding of astrophysical backgrounds, 
DD experiments can ultimately outperform the LHC for these mediator types (albeit with several model dependent 
caveats which render the results complementary to one anther in many regards). On the other hand, for the FCC 
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FIG. 8: 100 TeV exclusion limits on mediator mass versus dark matter mass for a given coupling g based on multijet measure¬ 
ments. We show results for vector (upper left panel), axial-vector (upper right panel), scalar (lower left panel) and pseudoscalar 
(lower right panel). Vector and axial-vector mediators plots indicate the regions where the minimal mediator widths do not 
exceed mMEo/S, thus reducing the available parameter space reach to ttidm uimed/S and cutting off large values of the 
couplings. For (pseudo)-scalar mediators the width we use is always sufficiently narrow. 


large regions of parameter space can be accessed for the the axial and scalar mediator which cannot be accessed by 
DD experiments limited by the neutrino wall. This extension of the parameter space for these types of models is a 
significant benefit of the FCC since it is unlikely this parameter region will ever be reached by current generations of 
experiments. Finally we note that as the couplings grow larger the width of the mediator grows, eventually for large 
enough values of the coupling it can exceed tomed/2 and a particle interpretation of the mediator is no longer valid. 
We highlight these regions for the vector and axial vector case. For the scalar case, which corresponds to our multi-leg 
on-shell sample the limit is trivially imposed by the on-shell condition (for the region where todm > ’t^med/ 2). We 
therefore only present the contours for the cases where they are relevant. 

In Figures]^ and [T^ we show these exclusion contours for the characteristic to a simplified model analysis fixed value 
of the mediator couplings, (7 dm = Ssm = 1 for all 4 mediator models of Eqs. 0-0- To enable the direct comparison 
between different experiments/techniques, these figures show all five exclusion contours - the 14 TeV and the 100 TeV 
limits, using both the one-jet and the multi-jet analysis, together with the DD/ID non-collider limits/projections. 

It is interesting to note the dependence of the DD limits in the scalar mediator case on the number of quark degrees 
of freedom it couples to. Unlike the production mechanism at collider searches which is sensitive only to the heavy 
top quark, the DD limits are sensitive also to light degrees of freedom thanks to the cancellation of the quark mass 
in the yq/iriq factor in Eq. (12). Thus, the DD limits are quite sensitive to choice of flavors that mediator couples to 
in the simplified model. The magenta contour in in Figure represents the inclusion of interactions with all quark 
flavors (as in the simplified model in Eq. 0). For a different choice of the simplified model, for example with only the 
top quark couplings to the mediator, the DD contour is shown in red. The difference between the red and magenta 
contours in the scalar mediator case in Fig. shows the sensitivity of the DD limits to a range of simplihed models; 
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FIG. 9: Mass limits for vector mediator models (left panel) and axial-vector models (right pannel) at 14 and 100 TeV colliders 
using the multi-leg and a single-leg analysis. We also show the neutrino wall limit of the direct detection. 




FIG. 10: Mass limits for scalar mediator models (left panel) and pseudo-scalar models (right pannel) at 14 and 100 TeV colliders 
using the multi-leg and a single-leg analysis. The neutrino wall affecting the direct detection experiments is shown in the left 
plot and the indirect detection limit for pseudo-scalars using FERMI-LAT data is shown as a tiny speck in the lower left 
of the plot on the right. 


at the same time the collider searches are are primarily sensitive to the scalar-to-top couplings^. For this parameter 
choice we note that the collider constraints lie below the neutrino wall for 1 ab“^, as the FCC collects more data the 
wall can be breached. As an example we plot the expected limit given 100 ab“^ of FCC data for the scalar mediator. 

Finally in Figure 11 we show the plots in terms of the spin-dependent and the spin-independent DM-neutron cross 
sections for a more traditional comparison of collider limits in terms of our simplified models with the limits/projections 
from the direct and indirect detection experiments. We compare the results in the tr, mjjM plane. The general pattern 
of Figs. 9p0 is reproduced, with the i^-wall for the spin-independent cases providing the strongest projected bounds. 
For the axial-vector and scalar mediators(with TeV-scale mediator masses as chosen in Fig. 11) our LHC contours 
cross the neutrino wall limits of direct detection experiments for todm ^ 10 GeV. At 100 TeV we hnd that collider 
bounds for the axial mediator are the strongest and universally below the z/-wall limit of direct detection, whilst for 
these parameter choices the scalar mediator and the j^-wall are comparable. In the pseudo-scalar case the last plot in 


® We note that in the previous figures the i/-wall curve corresponds to the magenta curve. 
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FIG. 11: Collider exclusion contours interpreted in terms of the spin-dependent and spin-independent cross sections plotted 
as functions of the dark matter mass. For the vector, axial-vector and scalar mediator models indicate the neutrino wall 
limit of the direct detection experiments. For the pseudo-scalar mediator model we show the indirect detection limits (using 
FERMI-LATdata [68]). 


Fig. [^demonstrates that both 14 TeV and 100 TeV collider bounds provide a multi-order of magnitude improvement 
over the current ID reach. 


IV. CONCLUSIONS 

We have presented a comprehensive study of the forthcoming and future hadron collider limits and projections at 14 
and 100 TeV for searchers of new physics associated with the Dark Matter sector. The dark sectors are characterised 
in this work in terms of four generic classes of simplified models Eqs. 0-0 where interactions between the Standard 
Model partons and the ‘invisible’ dark matter sector particles are described by four basic types of mediators: scalar, 
pseudo-scalar, vector or axial-vector particles. The dark matter particles we consider are produced via s-channel 
mediator exchange, see Fig. For collider searches of dark particles to be effective, two body decays of mediators 
produced on-shell should be kinematically possible, which requires that tomed > 2mDM, as can be seen in e.g. Figs.j^ 
It is then the ability to produce the mediator particle itself that underlines the efficiency of collider searches for dark 
matter, and not so much the particular species of Dark Matter the mediator decays into, in so far as these decays 
are kinematically allowed. Importantly, there is no requirement that once produced, the mediators should decay 
predominantly into the cosmologically stable Dark Matter, instead (and arguably more plausibly) they can have 
significant branching ratios for decay into any dark sector particles which are long-lived on collider scales. Hence we 
do not impose the relic density constraints on the dark particle production in our simplihed model treatment. 

Collider limits on the signal cross-sections for the cases of the 14 TeV LHC and the 100 TeV FCC are summarised 
in Fig. j^for all four mediator types. From these we have determined the collider reach on the (moM, wmed) mass 
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plane at 14 TeV in Fig.j^and 100 TeV in Fig[^ Both these sets of figures compare the collider reach with the neutrino 
wall limit of direct detection experiments and the current data from the indirect detection (in the latter case for the 
pseudo-scalar mediator models). We conclude that for scalar and axial-vector mediators collider searches are highly 
competitive and also complimentary to direct detection experiments, while in the pseudo-scalar case the collider limits 
are unchallenged. These figures also show the required strength of the DM coupling of the mediators, which in all 
cases remains largely in the perturbative regime. 

Figures and [T0| combine the 14 TeV and 100 TeV projections for collider searches and compare them with DD 
and ID limits at fixed values of mediator couplings 5dm = 5sm = 1- For vector and axial-vector mediators mediator 
masses of up to ^ 3 TeV can be probed at the LHC and extended to nearly 13 TeV at the 100 TeV FCC. For scalar 
and pseudo-scalar mediators, the 14 TeV LHC reach is of ttimed ~ 700 GeV and up to ~ 4-5 TeV for the FCC. Thus 
the Future hadron Circular Collider would be able to truly probe the few-TeV scales of the Dark Matter sectors, and 
in this regime it crosses over and goes beyond the projected limits of DD and ID experiments as can be seen from 
Figs. 9|T0 on the moM, Mmed mass plane and the plots in Fig. 11 which interpret collider searches in the language of 
DD experiments in terms of the spin-independent and spin-dependent cross-sections. Our FCC limits were presented 
assuming a rather modest data set of 1 ab“^, over the lifetime of the machine a much larger data set should be 
collected, and the limits/discovery sensitivity will be significantly enhanced. 

In summary, although our results for the FCC are somewhat speculative, it is clear that there is a huge potential for 
such a machine to probe a large parameter space for a variety of dark sector scenarios. Indeed an FCC style machine 
can provide access to regions of parameter space for axial and scalar mediators which lie beyond the potential limits 
of the cosmic neutrino wall. Such a possibility is an exciting aspect of the emerging physics program of the FCC, and 
in our opinion represents a strong motivation to press ahead with its ultimate construction. 
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